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^ ■ ABSTRACT 

Context. Modeling of pre-main sequence (PMS) stars through asteroseismology of PMS p-mode pulsators has only recently become 
possible, and spacebased photometry is one of the important sources of data for these efforts. We present precise photometry of the 
pulsating Herbig Ae star HD 142666 obtained in two consecutive years with the MOST (Microvariability & Oscilations of STars) 
satellite. 

Aims. Previously, only a single pulsation period was known for HD 142666. The MOST photometry reveals that HD 142666 is multi- 
periodic. However, the unique identification of pulsation frequencies is complicated by the presence of irregular variability caused by 
the star's circumstellar dust disk. The two light curves obtained with MOST in 2006 and 2007 provided data of unprecedented quality 
to study the pulsations in HD 142666 and also to monitor the circumstellar variability. 

Methods. Frequency analysis was performed using the routine SigSpec and the results from the 2006 and 2007 campaigns were then 
CNJ ' compared to each other with the software Cinderella to identify frequencies common to both light curves. The correlated frequencies 

were then submitted to an asteroseismic analysis. 
OO , Results. We attribute 12 frequencies to pulsation. Model fits to the three frequencies with the highest amplitudes lie well outside the 

. uncertainty box for the star's position in the HR diagram based on published values. Some of the frequencies appear to be rotationally 

split modes. 

Conclusions. The models suggest that either (1) the published estimate of the luminosity of HD 142666, based on a relation between 
, circumstellar disk radius and stellar luminosity, is too high and/or (2) additional physics such as mass accretion may be needed in our 

,] ■ models to accurately fit both the observed frequencies and HD 142666's position in the HR diagram. 

2 , 

Key words. Stars: pre-main sequence, (Stars:variables:) S Set, Stars: individual :HD 142666, Techniques: photometric 



1. Introduction 

A star evolving from the birthline to the zero-age main sequence 
(ZAMS) derives most of its energy (half of which heats the star 
and half of which radiates away) from the release of gravitational 
potential energy as the star collapses. These pre-main sequence 
(PMS) objects are characterized by observational features typ- 
ical for their early evolutionary stage, such as emission lines, 
infrared and/or ultraviolet excesses, (ir)regular brightness varia- 
tions, etc. The resulting variability of PMS stars is observed on 
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a wide range of time scales and amplitudes. Variations on time 
scales of weeks with amplitudes of a magnitude and more orig- 
inate from the circumstellar environment that contains gas and 
dust remnants of the stellar birth cloud. Accretion or chromo- 
spheric activity induces variations on time scales from several 
hours to days. Oscillations with periods of few hours to half an 
hour and amplitudes at the millimagnitude level are due to 6- 
Scuti-like pulsations in those PMS stars with the right combina- 
tion of mass, temperature and luminosity. 

There are two distinct classes of PMS objects: T Tauri stars 
and Herbig Ae/Be stars. The low-mass (< IMq) T Tauri stars 
are too cool (spectral types from late F to M) to be unstable to 
(5-Scuti-type p-mode pulsation. The Herbig Be stars are too mas- 
sive (4 to 10 Mq) for an observable PMS phase: it is believed that 
by the time they become first visible optically, they are already 
burning hydrogen in their cores. The intermediate-mass (1.5 to 
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4 Mq) Herbig Ae stars - on the other hand - cross the instabil- 
ity region in the Hertzsprung Russell (HR) diagram on their way 
to the main sequence. They have the right combination of mass, 
temperature and luminosity to be pulsationally unstable. 

Several pulsating PMS stars have been discovered and ana- 
lyzed within the last few years. A summary of their properties is 
given by Zwintz (2008 ). The periods of PMS pulsators are in the 
domain of classical 6 Scuti stars, i.e., between 20 minutes and 6 
hours, and their amplitudes lie at the millimagnitude level. Pre- 
and post-main sequence stars of the same mass, effective temper- 
ature and luminosity mostly differ in their inner structures, but 
their atmosphere properties are quite similar (Marconi & Palla 
11998 ). PMS stars lack regions of already processed nuclear ma- 
terial and seem not to be affected by strong rotation gradients 
which can complicate their inner structures. 

Asteroseismology is the only method to probe the interiors 
of these stars, thanks to the fact that the excited pulsational in- 
stabilities depend strongly on the stellar density profile (Suran 
et al. l2001l l. Since the higher-order frequency spacings for pre- 
and post-main sequence stars are different (Suran et al. 120011 1, 
the evolutionary phase of a field star may be determined from its 
pulsational eigenspectrum. The first successful seismic models 
of PMS stars have already been obtained (e.g., Guenther et al. 
120071 Zwintz et al. 2007 ). 

In this paper we discuss MOST spacebased observations of 
the pulsating Herbig Ae star HD 142666 and present a first as- 
teroseismic model for this star. 



2. HD 142666 

The intermediate-mass Herbig Ae star HD 142666 (V 1026 Sco) 
has been studied frequently in the past. Its irregular brightness 
variations, accompanied by color changes (i.e., the star becomes 
redder when fainter), can be explained by a circumstellar dust 
disk seen nearly edge-on (e.g., Meeus et al. 1998 ) with grains 
smaller than -1 fim (Lecavelier des Etangs et al. 2005 ). When 
the dust clouds are in the observer's line-of-sight, they obscure 
the star and cause the reddening. This phenomenon of 'non- 
periodic Algol-like brightness minima' (also called "UX Ori 
type variations") is typical for Herbig Ae/Be stars (Grinin et al. 
|1994t . Even at maximum brightness HD 142666 shows a signif- 
icant amount of reddening (Meeus et al. |1998l l, which indicates 
that the light from the star is always diluted by its dusty circum- 
stellar environment. Note that this also makes difficult a reliable 
determination of the absolute magnitude and, hence, the lumi- 
nosity of the star. 

HD 142666 is a Type II Herbig Ae star; i.e., it shows a 
double-peaked H„ emission profile (e.g. Vieira et al. |2003l l. Its 
projected rotational velocity, v ■ sin / , was determined to lie be- 
tween (Meeus et al. 119981 1 and 91 kms^^ (Vieira et al. 
|2003ll. With a spectral type of A8 Ve and log = 3.857 (Vieira 
et al. 120031 ), HD 142666 is perfectly suited to search for, and in- 
vestigate, PMS pulsation. 

A single pulsation frequency of 2 1 .43 d^"^ (i.e., period of 1 . 12 
hours) was discovered by Kurtz & Miiller ( 120011 ) based on 6.6 
hours of observations obtained during a single night. They iden- 
tify this frequency as a fundamental mode according to the cal- 
culation of the pulsation constant, Q. For a detailed pulsational 
analysis and asteroseismic study, longer time series of data with 
high time sampling and excellent photometric precision were es- 
sential. Hence, HD 142666 was selected as a MOST Primary 
Science Target. 



3. MOST observations 

The MOST space telescope (Walker et al. 120031) was launched 
on 30 June 2003 into a polar Sun-synchronous circular orbit of 
altitude 820 km. (MOST's orbital period is 101.413 minutes, 
corresponding to an orbital frequency of ~14.2 t/"'.) From its 
orbital vantage point, MOST can obtain uninterrupted observa- 
tions of stars located in its Continuous Viewing Zone (CVZ) 
for up to 8 weeks. The MOST satellite houses a 15-cm Rumak- 
Maksutov telescope feeding a CCD photometer through a single, 
custom broadband optical filter (covering wavelengths from 350 
to 750 nm). 

MOST can supply up to three types of photometric data si- 
multaneously for multiple targets in its field. The mission was 
originally intended only Fabry Imaging, in which an in-focus 
image of the entrace pupil of the telescope - illuminated by a 
bright target star (V < 6) - is projected onto the instrument's 
Science CCD by a Fabry microlens (see Reegen et al. 2006 for 
details). After MOST was operating in orbit, the pointing perfor- 
mance of the satellite was improved so much that a new mode of 
observing. Direct Imaging, was made practical. Direct Imaging 
is much like conventional CCD photometry, in which photom- 
etry is obtained from defocussed images of stars in the open 
area of the CCD not covered by the Fabry microlens array field 
stop mask. In the original mission, no scientific information was 
available from the guide stars used for the ACS (Attitude Control 
System), but now precise photometry is possible for these stars 
as well (see, e.g.. Walker et al. 2005 and Aerts et al. 2006). 

Due to its brightness of V = 8.81 mag, the PMS pulsator 
HD 142666 is best suited for the Direct Imaging mode of pho- 
tometry. The star lies slightly outside the MOST CVZ, so it can 
be observed only for a part of each 101 -min orbit. But even 
this time coverage represents an extremely high duty cycle com- 
pared to what is obtained with groundbased measurements. On- 
board exposures are 1.5 seconds long (to satisfy the cadence of 
guide star ACS operations), but 14 consecutive exposures are 
"stacked" on board to produce integrations 21 sec long, sampled 
about 3 times per minute. 

MOST observed HD 142666 for the first time during 22 
April - 4 May 2006, for 11.5 days (see Figure[T]top panel) with a 
duty cycle of 34%. (Note that the duty cycle represents the frac- 
tion of each 101 -min orbit covered by the data; there are no long 
gaps in the time series.) The irregular UX Ori-type variations 
with amplitudes of nearly one magnitude due to the circumstel- 
lar dust disk are clearly visible in the light curve. 

HD 142666 was reobserved by MOST during 18 April - 28 
May 2007 (see Figure[T]middle panel) with a duty cycle of 39%. 
The total time span of the 2007 observations (~39 days) is more 
than 3 times that of the 2006 data. The gap in the light curve 
between t - 2696 d and t - 2697 d results from an interruption 
of the HD 142666 observations for a high-priority MOST target 
of opportunity. The bottom panel in Figure [1] is an enlargement 
of the 2007 light curve where the pulsational variations at the 
millimagnitude level can be seen clearly. 

4. Data Reduction & Frequency Analyses 

For reduction of MOST Direct Imaging photometry, two dif- 
ferent methods were developed. The first, by Rowe et al. 
(I2006al I2006bl ). combines classical aperture photometry and 
point-spread function fitting to the Direct Imaging subrasters. An 
independent data reduction pipeline for spacebased open field 
photometry, which also includes automated corrections for cos- 
mic ray hits and a stepwise pixel-to-pixel decorrelation of stray 
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Fig. 1. MOST time series photometry of HD 142666: 2006 observations (top panel) and 2007 observations (middle panel) on the 
same scale. A zoom into the 2007 light curve shows the pulsational variability (bottom panel). 



light effects on the CCD, was developed by Huber & Reegen 
(120081 1. Both routines were applied to both sets of HD 142666 
photometry, and no significant differences in the qualities of the 
reduced light curves could be identified. In the following anal- 
yses, the plots show data reduced using the Rowe et al. (I2006al 
12006b) method. 

The main challenge in the frequency analysis of HD 142666 
was the careful distinction of the large-amplitude variations 
caused by the circumstellar material from the small-amplitude 
pulsations. SigSpec (Reegen [20071 1 and Cinderella (Reegen et 
al.|2008) were used in combination to identify the pulsation fre- 
quencies. 

SigSpec (Reegen r2007l l computes significance levels for am- 
plitude spectra of time series with arbitrary time sampling. The 
probability density function of a given amplitude level is solved 
analytically and the solution includes dependences on the fre- 
quency and phase of the signal. 

Cinderella (Reegen et al. I2008I I is based on the principles of 
SigSpec and compares sets of frequencies from a target and one 
or more comparison light curves. A frequency is considered to be 
present in all data sets (and, hence, to be instrumental) if a target 
peak appears at the same frequency (within the frequency reso- 
lution) with at least the same significance as in the comparison 
data set(s). Such a coincident peak is then assigned a negative 
conditional significance. 



For each of the two HD 142666 light curves and their cor- 
responding background light curves, amplitude and significance 
spectra were computed. The results of these calculations were 
then submitted to Ce^jderella separately for each year of the 
observations, to identify coinciding peaks between the target and 
the background. Such peaks can originate from the background 
itself (due to scattered Earthshine variations) or from instrumen- 
tal effects. The irregular light variations of HD 142666 intro- 
duce a 'pseudo-periodicity' and an excess of significant peaks 
in the low-frequency domain (i.e., from to 3 li ') which can 
be seen, e.g., in Figure |4] The frequency range where 5-Scuti- 
like pulsations are expected - i.e., from ~5 to 70 t/"' - is well 
separated from the peaks caused by the irregular variability. 
Nevertheless, the formally significant frequencies in this domain 
have to be interpreted carefully for several reasons. 

First, alias frequencies from the low-frequency domain can 
appear at higher frequencies due to the spectral window. A man- 
ual check against such spurious peaks as well as a dedicated test 
(see Section |4!4] i were carried out. Second, instrumental frequen- 
cies related to the orbit of the satellite, its harmonics and Ic/"' 
sidelobes are present in the range of the expected pulsation fre- 
quencies. All significant peaks that can be related to one of the 
above mentioned instrumental frequencies within the frequency 
resolution (computed according to Kallinger et al. i2008i were 
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Fig. 2. Frequency analyses of the 2006 MOST photometry of HD 142666. The amplitude spectrum is shown at top left (labelled "a") 
and the spectral window at top right. Panels "b" to "f show power spectra where the identified pulsation frequencies are marked 
with arrows. Panel "g" (bottom right) is the residual power spectrum after prewhitening of all 151 formally significant frequencies in 
the data. The solid grey lines mark the MOST satellite orbital frequency and its harmonics. The dotted grey lines are Ic/ ' sidelobes 
of these frequencies. Note the different scalings on the y-axes. 



discarded. Finally, the identified pulsation frequencies need to 
be significant in the analyses of both HD 142666 time series. 

4.1. 2006 data 

The top left panel (labelled "a") in Figure|2]shows the amplitude 
spectrum of the 2006 data from to 100 li"' where the solid grey 
lines mark the orbital frequency of MOST and its multiples and 
the dotted grey lines represent the corresponding It/ ' sidelobes. 
The spectral window is shown in the top right panel. 

We find 151 peaks (grey lines in the top layer of Figure 
m to be formally significant; i.e., their significance criterion is 
larger than 5 (Reegen|2007]l which corresponds to an amplitude 
signal-to-noise (S/N) ratio > 4 (Breger et al. 119931 Kuschnig et 
al. 1997). About half of these frequencies lie between and 3 
c/"', and can be attributed to the irregular variability caused by 
the dusty circumstellar environment. Another large fraction of 
the formally significant peaks can be related to the orbital fre- 
quency of the satellite, its harmonics and Ic/"' sidelobes. 

In panel "b" of Figure[2l the appearance of alias frequencies 
due to the gaps in the data per MOST orbit can be seen: aliases of 
the two highest pulsation frequencies at 21.250 d^^ and 22.009 
d^^ appear at +I4.2d^^ which is the orbital frequency of MOST. 
After prewhitening these two main frequencies, the respective 
aliases also disappear. 



The residual noise level in the 2006 data set is 38 ppm. 



4.2. 2007 data 

The top left panel (labelled "a") in Figure[3]shows the amplitude 
specti-um of the 2007 data from to 100 c/ '. Since the 2007 
MOST time series is more than 3 times longer than the 2006 
data set, the amplitude spectrum of the 2007 photometry has an 
obviously lower noise level. 

The regularly spaced peaks between and 14 c/"' visible in 
the amplitude spectrum are also present in the spectral window 
(top right in Figure[3]). The cause is gaps in the data due to inter- 
ruptions to observe yet another high-priority MOST target of op- 
portunity, monitored every half-day during the 2007 HD 142666 
run. This regular spacing of 2 d^^ and its harmonics appear in 
the amplitude spectrum of the background light curve with even 
higher amplitudes. When submitting the star and the background 
light curves to Cinderella, these frequencies get negative con- 
ditional significances and are therefore rejected. Figure |5] shows 
the results of the analysis with Cinderella. 

There are 451 formally significant frequencies in the 2007 
MOST data (grey lines in the bottom panel of Figure|4]i, where a 
third of the peaks lie between and 3 li ' and are due to the irreg- 
ular light variations of the star. As with the 2006 data, the alias 
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Fig. 3. Frequency analysis of the 2007 MOST observations of HD 142666. The format of the figure is the same as Figure 1. Panel 
"g" (bottom right) is the residual spectrum after prewhitening of all 415 formally significant frequencies found in this time series. 
Note the different scalings on the y-axes. 



frequencies due to the spectral window can be seen at ±\A.2d ' 
around 21.252 d^^ and 22.008 t/"' (see Figure[3l panel "b"). 
The residual noise level in the 2007 data set is 30 ppm. 

4.3. Pulsational frequencies 

From the numerous formally significant peaks in the amplitude 
spectra of HD 142666, 12 were identified to likely originate from 
pulsation because they meet the following criteria: (a) They are 
not related to the orbital frequency of the satellite, its harmonics 
or l(i ' sidelobes. (b) They appear significantly in the data sets 
from both years, (c) They cannot be attributed to the peaks intro- 
duced by the irregular light variations and their alias frequencies, 
(d) They passed the test described below in 14.41 

The 12 frequencies, their amplitudes and significances are 
Usted in Table [1] where the errors given in parentheses are 
computed using the formulae given by Kallinger et al. ( 120081 ). 
Although the significances and corresponding amplitudes of 
some of the listed frequencies are quite low, they all meet the 
criteria described above. 



4.4. Test with 2007 subset iiglit curve 

The reliability of the selected 12 pulsation frequencies was in- 
vestigated in the following way. The longest part of the light 
curves with minimal irregular light variations lies in the 2007 



data set from t - 2668.9 li to f = 2613.9 d (see top panel in 
Figure |6ll. 

First, a frequency analysis of this 5-day subset light curve 
was carried out with SigSpec. Although this is the "quietest" part 
of the light curve, irregular variability is still clearly present. To 
reproduce the shape of the irregular light variations in this data 
subset, moving averages were computed. Using a time interval 
of 4 hours for the boxcar allows one to smear out the pulsation 
effects deliberately, without significant suppression of the irreg- 
ular variations. The bottom panel of Figure |6] shows the subset 
light curve (grey dots) where the boxcar is overplotted (black 
dots). 

A frequency analysis of the moving average was then com- 
puted with SigSpec. The results of the subset light curve and 
the boxcar light curve were then compared using Cinderella. 
If frequencies are found in both data sets (i.e., they have nega- 
tive conditional significances in the output of Cinderella), they 
very likely originate from the circumstellar environment and are 
not caused by pulsation. 

All 12 frequencies previously attributed to pulsation were 
only present in the computation of the subset light curve. They 
do not show up in the frequency analysis of the boxcar and have 
positive conditional significances in Cinderella (see Figure|7]i. 
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Table 1. Pulsation frequencies of HD 142666 identified in both data sets (2006 and 2007), their amplitudes and significances, where 
the respective frequency errors are given in parentheses. The table is sorted by increasing frequency, where the frequencies are 
numbered according to their significances. 
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Fig. 4. Significance spectra of the 2006 data (top) and 2007 data 
(bottom). Grey lines mark all formally significant peaks; black 
lines are the 12 identified pulsational frequencies common to 
both data sets. 



5. Asteroseismic Analysis 

5.1. Model grid 

We compare the observed frequencies from the 2007 run for 
HD 142666 to asteroseismic models in order to determine the 
nature of the pulsations and to constrain the star's effective tem- 
perature, luminosity and mass independently. We note that the 
model fits to the 2006 frequencies (and averaged 2006 and 2007 
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Fig. 5. A closer look at the low-frequency domain of the results 
of the 2007 data analysis with Cinderella. The 2-(i ' frequency 
is known to be due to the gaps introduced in the HD 142666 
time series. This and other peaks with negative values of their 
conditional significance are shown by grey lines. Peaks that ap- 
pear only in the stellar light curve have positive conditional sig- 
nificance (black lines). 



frequencies) are nearly identical, which is as expected since the 
differences between corresponding frequencies is less than 1.0 
fjHz, and typically less than 0.25 yuHz. 

Grids of models were constructed using the Yale Rotating 
stellar Evolution Code (YREC; Guenther et al. T^92]l. P MS 
models were evolved from the Hayashi track (Havashi |l961| l be- 
fore deuterium burning to the ZAMS. The grids include models 
with masses ranging from 1.00 to 5.00 in steps of 0.01 Mq. 
Each evolutionary track is resolved into approximately 1000 
models. A near-solar composition (Z = 0.02, Y = 0.27) was as- 
sumed. 

The constitutive physics of the models are current and in- 
clude OP AL98 (Iglesias & Rogers .1996 ), the Alexander & 
Ferguson (|1994| l opacity tables, and the Lawrence Livermore 
National Laboratory equation of state tables (Rogers 119861 
Rogers et al. 1996). The mixing length parameter used to de- 
scribe the temperature gradient in convective regions according 
to the Bohm-Vitense (il958) mixing-length theory, was adjusted 
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Fig. 6. Five-day subset light curve from the 2007 data set (top) 
and a moving average fit over a time interval of 4 hours (in black) 
plotted over the light curve (bottom). 
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Fig. 7. Comparison between the frequencies of the 5-day sub- 
set light curve and the moving average fit over a time interval 
of 4 hours using Cinderella. The grey lines show the output of 
Cinderella, where peaks with negative conditional significances 
appear in both data sets, and are unlikely to be caused by pulsa- 
tion. The 12 pulsation frequencies identified in HD 142666 have 
positive conditional significances. Their positions are shown as 
black lines for better visibility. 



from calibrated solar models. The PMS evolutionary tracks used 
in this analysis compare well with the tracks of D'Antona & 
Mazzitelli ( |1994| l. Our models do not include mass accretion 
(Palla & Stahler T"992| 11993 '), which may occur prior to initial 
deuterium burning. We used Guenther's nonradial stellar pulsa- 
tion code (Guenther |1994| ) to compute the adiabatic pulsation 
spectra. 

We used a dense and extensive grid of model spectra to find 
the best match to the observed frequency spectrum. This method 
was originally developed by Guenther & Brown ( I2004I I for mod- 
eling stars in more advanced evolutionary stages where mode 
bumping complicates the search for the best match to the ob- 
served oscillations. The quality of the match is quantified by the 
simple relation as described, e.g., in Guenther et al. ( 120071 ), 
that compares the model and observed frequencies weighted by 
the model and observational uncertainties. The grid of oscilla- 
tion spectra is searched to locate local minima in chi-squared. A 



value of < 1 means that the difference between the observed 
frequencies and the corresponding nearest model frequency is 
less than the rms of the model and observational uncertainties. 
We assumed the observed frequencies have an uncertainty of 
+0.3juHz (corresponding to the reciprocal of the 2007 observ- 
ing run duration). 



5.2. Modelling HD 142666 

Before trying to find a model whose oscillation spectrum 
matches an observed set of frequencies, we first preview the 
observed frequencies in an echelle diagram. This enables us 
to quickly check if the observed frequencies show the ridge- 
like structure expected for p-modes. The frequencies of p- 
modes, especially those of PMS stars whose chemical composi- 
tions are nearly homogeneous, are relatively uniformly spaced. 
Specifically, we expect modes of similar /-value to be sepa- 
rated by the large spacing, defined as the frequency difference 
between n and n + 1 p-modes (of similar Z-value.) By plotting 
the frequency of a mode versus its frequency modulo the large 
spacing in an echelle diagram, one can immediately identify se- 
quences of frequencies or ridges that correspond to similar l- 
valued modes. The large spacing is approximately equal to the 
characteristic frequency spacing (Tassoul 1980) which can be 
computed directly from a stellar model of the star. 

We found that no single folding frequency for the echelle di- 
agram allows all 12 observed frequencies to fall onto /-ridges. 
Indeed, it was immediately clear, based on their proximity, that 
the two frequencies f7, f8 and the three frequencies f5, f9, and 
fl2 appear to be rotationally split modes. HD 142666 has a pro- 
jected rotational velocity v ■ sin / of about 50-100 A:m5 ', which 
is consistent with this interpretation. Because our model fre- 
quencies are obtained from a non-rotating stellar model, no ro- 
tationally split mode frequencies are included in our grid search. 
Hence, the oscillation spectrum searching algorithm will fail to 
fit, or be skewed by, any observed frequency that is an m + 
(i.e., not the central) split mode. Consequently, we included only 
the central frequency, f7, and f9. The rotationally split assump- 
tion is, of course, speculative. Regardless of their true nature, as 
can be seen in the echelle diagram (FigurefTTT). we are unable to 
find any model fits to all these frequencies simultaneously since 
their spacings does not match simultaneously the small and large 
spacings, as would be required to fit all three frequencies. 

Although the ultimate goal of the seismic analysis is to de- 
termine a model of HD 142666 based on observed frequen- 
cies alone, it is extremely useful to have some prior information 
about the star. Unfortunately, because HD 142666 is surrounded 
by a circumstellar disk, its location in the HR-diagram is uncer- 
tain. The effective temperature of HD 142666 as found in the lit- 
erature varies wildly. For example, Vieira et al. (I2003I I find T^s 
= 7200 K while Dominik et al. (2003 ) find T^s = 8500 K for this 
star. The recent analysis of data obtained with the ESPADONS 
spectrograph (J. Grunhunt & G. Wade, priv. comm.) yielded red 
- 7500 K. Here we adopt for reference the value of Vieria et al. 
(120031 ), log Teff = 3.857±0.05. The only estimate of luminosity 
we could find in the literature is by Monnier et al. (2005) who 
obtain L - 8.8 ± 2.5 derived from a relation between cir- 
cumstellar disk radius and stellar luminosity. The result is highly 
uncertain, since the input parameters of the disk are, themselves, 
uncertain. 

We computed the characteristic frequency for the models in 
our PMS grid that are in the vicinity of HD 142666's location in 
the HR diagram and plotted them in the contour plot shown in 
Figure [8] Figure [8] shows that if the star is correctly positioned 
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Fig. 8. Contour plot of the characteristic frequency spacing (in Fig. 9. Contour plot of the frequency (in //Hz) of the fundamental 
IdUz) near the star's position in the HR-diagram showing that radial mode near the star's position in the HR-diagram. 
the large spacing is expected to be between 45 and 80 fjHz. 



in the HR diagram, then we can expect the large spacing to be 
between 45 and 80 //Hz. As we show below, our best model fits 
have a much smaller large spacing (~15 //Hz). 

Figure|9]is a contour plot on which we show the frequency of 
the fundamental radial model for all the models in our PMS grid 
near HD 142666's presumed location in the HR diagram. If the 
HR diagram location of HD 142666 is correct, then the p-modes 
must have frequencies greater than -100 //Hz. Consequently, 
without further analysis we can immediately see that, if the ob- 
served frequencies down to 67 //Hz are all intrinsic p-modes, 
then HD 142666's presumed luminosity and/or temperature are 
incorrect. Of course, it is also possible that, if intrinsic, the low- 
est frequencies are g-modes. Without corroborating oscillation 
observations, and with a highly uncertain HR diagram location, 
we cannot say for certain. 

We began our spectrum fitting analysis by considering only 
the three most significant frequencies fl, f2, and f3. In the HR- 
diagram shown in Figure [TO] we plot all the models that have 
oscillation frequencies that match the three observed frequen- 
cies within the assumed frequency uncertainty of +0.3 //Hz. The 
model oscillation spectra include / = 0, 1, and 2 p-modes. We 
find a variety of models spread throughout the HR-diagram that 
fit the three most significant frequencies with < 1 . None, 
though, lie near the HR diagram position of the star. 

No model spectra in our PMS grid were found to fit, with 

< 1, all the observed frequencies excluding f5, f8, and fl2 
(the frequencies tentatively identified as side lobes to rotation- 
ally split modes). We, therefore, considered the possibility that 
some of the frequencies are not p-modes. After investigating 
several combinations of frequencies, we applied the following 
restrictions to obtain the final five frequencies that we are able 
to fit. First, we considered only frequencies > 125 //Hz, thus 
eliminating frequencies f4, f6, flO and fl 1 from our fitting tests. 
Second, we took the central frequency only in the triplet f5, 
f9, fl2, i.e., f9. And, finally, we chose the most significant fre- 



quency, fl, from the doublet f7, f8. Figure [TOl shows the best 
model fits to the resultant set of frequencies fl, f2, f3, fl, and f9. 
The large filled circles correspond to model fits with;if^ < 1 and 
the large open circles to models with < 2. The model fits to 
all five frequencies is a subset of the model fits to the three most 
significant frequencies. 

Figure [TT] is an echelle diagram of the model frequencies 
from the best model fit, i.e., the model with the lowest (~0.3), 
to the selected five frequencies. The model has log T^ff - 3.8267, 
log L/Lq- 1.6903, radius = 5.2Rq and mass = 2.96Mq . The five 
fitted frequencies are indicated by large filled circles. The re- 
maining observed frequencies are noted by large open circles. 
The / = model frequencies are shown with small open squares, 
the I - 1 mode frequencies are shown with small open triangles, 
and the / = 2 and 3 model frequencies are shown with small open 
circles running parallel to the I = and 1 modes, respectively. 
The echelle diagram shows that the five observed frequencies 
are well fit by the / = 1 and 2 p-modes. We note that we do not at 
any time during the fitting process force specific frequencies to 
be fitted, for example, by only / = modes. If we consider only 
models that lie within the estimated effective temperature range 
(ignoring the luminosity discrepancy) that fit the three most sig- 
nificant frequencies then we arrive at nearly the same models, 
and produce nearly the same echelle diagrams as our best model 
fit to all five frequencies. Essentially, the additional two frequen- 
cies help constrain the star's location in the HR diagram to coin- 
cide with the estimated effective temperature range. 

The identification of rotational splittings for the triplet f5, 
f9, and fl2 and the doublet f7 and f8 is somewhat uncertain in 
that the splitting spacings are not the same for the lower fre- 
quency pair compared to the higher frequency triplet. The larger 
splitting for the lower frequency pair (assuming that the missing 
frequency is not the central frequency) implies that the core of 
HD 142666 is rotating more rapidly than the surface layers. For 
our best fitting model a splitting on the order of 2 /<Hz, under 
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Fig. 10. The HR-diagram position of HD 142666 (filled square 
with error bars) together with all the models that have oscillation 
frequencies matching the combinations of fl, f2, and f3 (small 
filled circles) and fl, f2, f3, f7, and f9 (large open circles for 
< '1 and large filled circles ici'^ ^ 1) with the frequency 
uncertainty. PMS evolutionary tracks are also plotted to show 
the relative evolutionary phase of the model fits. The masses, in 
solar mass units, label the individual tracks. 



the assumption of solid body rotation, yields a surface rotation 
period of ~6 d and a surface rotation velocity at the equator of 
~50 kms^^. This is a typical rate for PMS stars and compares 
well to the observed v • sin / measurements within the uncertain- 
ties of our own velocity determination which are +50%, owing 
to the uncertainty in frequency splitting itself which is on the 
order of 1 jMz, i.e., 50%. 

If the five frequencies are p-modes and if our classically de- 
termined PMS models are reasonable representations of the inte- 
rior structure then the observed frequencies are able to constrain 
the mass, luminosity, and surface temperature of HD 142666. 
The asteroseismic mass, for example, is far away from the mass 
that would be determined from simply fitting the star's loca- 
tion in the HR-diagram, i.e., ~1.1Mq (see Figure [TOl). There 
are, though, several issues that need to be resolved. We need 
to understand why the best fitting models lie well away from the 
star's presumed location in the HR diagram. Possibly this could 
be due to the presence of the circumstellar disk, which makes 
photometric determinations of the star's luminosity and effec- 
tive temperature highly uncertain. We also have to wonder why 
none of the frequencies fit are radial modes (i.e., I = 0). Most 
likely, our mode identifications and the specific model fit are not 
completely correct. For example, if the outer envelope and at- 
mosphere of the models are slightly off, a reasonable possibility 
considering the amount of dust in the neighborhood of the star, 
then the higher-frequency modes that are sensitive to the outer 
layers of the model, will also be off. A model uncertainty of 
+ 1//Hz is possible, which is enough to allow / - modes to be 
fit, but not enough to push the HR-diagram position of the best 
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Fig. 11. Echelle diagram of the model from the best model fit 
to the five frequencies fl, f2, f3, 17, and f9. Fitted frequencies 
are plotted as large filled circles, frequencies not used in the fit 
as large open circles. I - modes are shown with small open 
squares, I = 1 modes as small open triangles, I = 2 and 3 modes 
with small open circles running parallel to the I - and 1 modes, 
respectively. 

fitted models toward the star's apparent position. We believe a 
reasonable next step is to include mass accretion, as is done in 
the best models of Palla & Stabler (1992, 1993), and see how 
large is the effect and if there is any improvement in the fits. 

At this time we note that adjustments to mixing length pa- 
rameter, the helium abundance, and heavy element abundance 
will perturb the location of the models in the HR diagram only 
by a small amount (e.g., ~300K) - far from enough to account 
for the gap between the star's putative location in the HR dia- 
gram and the seismically fitted models. 

6. Conclusions 

MOST high-precision time series photometry of the Herbig Ae 
star HD 142666 obtained in 2006 and 2007 was used to inves- 
tigate the pulsational properties of the star. This task was com- 
plicated by the fact that HD 142666 is surrounded by a dense 
circumstellar dust disk which causes irregular light variations in- 
troducing high-amplitude signal at low frequencies (i.e., from 
to 3 d^^). We were able to disentangle the pulsational variability 
with amplitudes at the millimagnitude level from the iiTegular 
variations with amplitudes of up to 200 mmag in the Fourier do- 
main. 

Although numerous peaks appear to be formally significant, 
we identify only 12 frequencies as pulsational in origin in the 
data sets of both years. Therefore, only these were used for the 
asteroseismic analysis. 

For our model fitting, we further restricted the set of frequen- 
cies by excluding frequencies below 125 fjHz. Also, based on 
their proximity in frequency space, we included only the central 
frequency from the f5, f9, fl2 triplet, and the most significant 



10 



K. Zwintz et al.: MOSJ photometry of the enigmatic PMS pulsator HD 142666 



frequency from the f7, f8 pair, since we suspect they are rotation- 
ally split modes. We, thus, restricted the fits to the five modes fl, 
f2, f3, f7, f9. 

We are able to obtain good fits (x^ < 1) to the frequencies for 
/ = 0, 1 and 2 p-modes. Significantly, though, all the model fits 
lie well outside the HR diagram uncertainty box for HD 142666. 
Even when we fit only the three most significant frequencies 
(fl, f2, and f3), the models with;\f^ < 1 lie far away from the 
presumed location of the star in the HR diagram. Furthermore, 
we note that the fundamental p-mode frequencies in the vicinity 
of HD 142666's location in the HR diagram are higher than the 
lowest frequencies observed. Higher luminosity models, though, 
do have lower fundamental frequencies. 

The failure to find a model that fits even the three most sig- 
nificant modes that at the same time lies within the uncertainty 
box for the star's location in the HR diagram leads us to suggest 
three areas to investigate further: 

1. The models are too crude and need additional modeling 
physics such as mass accretion during the pre-deuterium burn- 
ing phases to fit the frequencies and the star's location in the HR 
diagram simultaneously. The presence of the circumstellar disk 
supports this idea but the size of the effect on p-mode frequen- 
cies needs to be determined. 

2. The third mode in the possibly rotationally split group fl, f8 
needs to be identified to help support the hypothesis that the 
modes are rotationally split. This would also help mode iden- 
tification in that the modes would then have to be nonradial. 

3. The luminosity of the star needs to be better constrained. This 
would rule out either the model fits or the p-mode nature of the 
lowest frequencies observed. If the star's luminosity is correct 
then the lowest observed frequencies are either non-intrinsic or 
they are g-m odes. Note that for HD 142666 no HIPPARCOS 
(ESA 119971 ) parallax is available and the error of the parallax 
fr om the ASCC 2.5 Catalog (Kharchenko 2001 Kharchenko et 
al. 120051 ) is a factor of three larger than the parallax itself. 
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